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Abstract. The least action principle, through its variational formulation, possesses a final- 
ist aspect. It explicitly appears in the fractional calculus framework, where Euler-Lagrange 
equations obtained so far violate the causality principle. In order to clarify the relation be- 
tween those two principles, we firstly remark that the derivatives used to described causal 
physical phenomena are in fact left ones. This leads to a formal approach of irreversible 
dynamics, where forward and backward temporal evolutions are decoupled. This formalism 
is then integrated to the Lagrangian systems, through a particular embedding procedure. In 
this set-up, the application of the least action principle leads to distinguishing trajectories 
and variations dynamical status. More precisely, when trajectories and variations time ar- 
rows are opposed, we prove that the least action principle provides causal Euler-Lagrange 
equations, even in the fractional case. Furthermore, the embedding developped is coherent. 



1. Introduction 

The link between the least action principle and the causalilty principle has always been 
ambiguous. Poincare [17] sumed it up as follows: 

L'enonce du principe de moindre action a quelque chose de choquant pour 
I'esprit. Pour se rendre d'un point d un autre, une molecule materielle, sous- 
traite a Faction donnee de toute force, mais assujettie a se mouvoir sur une 
surface, prendra la ligne geodesique, c'est-d-dire le chemin le plus court. Cette 
molecule semble connaitre le point oil on vent la mener, prevoir le temps qu 'elle 
mettra a I'atteindre en suivant tel ou tel chemin, et choisir ensuite le chemin 
le plus convenable. L'enonce nous la presente pour ainsi dire comme un etre 
anime et libre. II est clair qu'il vaudrait mieux le remplacer par un enonce 
moins choquant, et oil, comme diraient les philosophes, les causes finales ne 
sembleraient pas se substituer aux causes efficientes^. 

Henri Poincare, La science et I'hypothese, 1902. 

Therefore, how is it possible to obtain causal equations, i.e. equations taking into account 
only the past states, by using a principle which depends on the whole temporal interval? 
What does the information on the future become? Moreover, several approaches have been 
developped [H HI EJ [20] to generalize the least action principle and the Euler-Lagrange equation 



^ The very enunciation of the principle of least action is objectionable. To move from one point to another, 
a material molecule, acted upon by no force, but compelled to move on a surface, will take as its path the 
geodesic line - i.e., the shortest path. This molecule seems to know the point to which we want to take it, 
to foresee the time that it will take it to reach it by such a path, and then to know how to choose the most 
convenient path. The enunciation of the principle presents it to us, so to speak, as a living and free entity. It 
is clear that it would be better to replace it by a less objectionable enunciation, one in which, as philosophers 
would say, final efi^ects do not seem to be substituted for acting causes |18) . 
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to the fractional case. In this formahsm, the derivatives are non-local, which makes the 
past and future of the functions appear explicitly. Because of the simultaneous presence 
of left and right derivatives, none of those equations respect the causality principle. This 
difficulty may have been seen as a definite failure, and alternatives have notably been proposed 
in [71 [22] to get around this problem. However, because we believe that the least action 
principle should remain fundamental in any formalism, we choose to keep this approach - 
particularly the one in [5] -, and we prove in this paper that causality may be respected. To 
this purpose, we formulate the following remark: when one observe a phenomena which one 
wants to describe using a differential equation, one only have access to the left derivatives of 
the functions, i.e. to the differential operators dependent on the past values of the function. 
Thus, this characteritic represents a trace of the time arrow, and the differential equation 
becomes attached to the forward temporal direction. If we assume the existence of a similar 
differential equation, but related to the backward temporal evolution, we obtain a formal 
approach of irreversibility. Difficulties about causality inherent to the least action principle 
can be solved using this formalism. More precisely, we show that using a new embedding, 
termed asymmetric embedding, it is possible to obtain causal Euler-Lagrange equations. In 
doing so, we observe that the information on the future lies in fact in the variations used by 
this variational method. The virtual status of these could hence moderate the finalist aspect 
of the least action principle. In addition, we prove that the asymmetric embedding is coherent, 
i.e. that this procedure is globally compatible with the least action principle. 

The formal approach on irreversibility is first adressed in section [21 The asymmetric embbe- 
ding is introduced in section [3] and applied to the Lagrangian systems in section [H This leads 
to a causal Euler-Lagrange equation, obtained in section [5l Application of this formalism to 
few examples of derivatives is reported in section [6] while those results are discussed in section 



2.1. Dynamics and causality. In physics, the causality principle means that the state of 
a system at a time t is completely determined by its past, i.e. by its states at times t' , 
t' < t. Therefore, if a system is described by variables x G M", the variations of those ones 
(for example the velocities) should only depend on the past instants. The derivatives, which 
express those variations, are hence left derivatives. They will be denoted by V~^. The following 
definition formalises these ideas. 

Definition 1. The evolution of a system is said causal in the direction past — > future if it 
can be written as 



m 



2. An approach of irreversibility 




(2.1) 



A first example of operator P"*" is the usual left derivative 



, , x(t) — x(t — e) 
d+x{t) = lim — ^ 

e^O ^ 



(2.2) 



e>0 



However, following [9], we postulate that the evolution of a general physical system is a 
priori irreversible: the dynamics in the direction future —>■ past cannot be described by (|2.1|) . 
In order to have a complete description of the system, a supplementary differential equation has 
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to be introduced which accounts for the evolution towards the past. The evolution operators 
are in this case right derivatives, denoted by P", and lead to the following definition. 

Definition 2. The evolution of a system is said causal in the direction future — > past if it 
can be written as 

/_ (x{t),V~x{t), {V-fx{t),t) = 0. (2.3) 
For the previous example, the operator is the usual right derivative 

d^t) = lim + 

Now we can precise our formal approach of irreversibility. 

Definition 3. A system is said reversible if (12. ip and (|2.3I) have the same solutions. Other- 
wise, it is said irreversible. 

Remark 1. We emphasize the formal aspect of this definition: our goal here is not to un- 
derstand the physical origin of irreversibility. This problem, related to Boltzmann's work on 
entropy in the 1870 decade, is still discussed in the physics community. We simply mention 
that recent answers have been proposed through chaotic systems |19l I23j . 

Even if our approach is formal, it will be useful to understand better the least action 
principle, more precisely its relation with causality. 

2.2. Asymmetric dynamical representation. In the rest of the paper, we will consider a 
system evolving in M", during a temporal interval [a, b]. This system will de denoted by <S. 
We introduce the vector space U defined by 

U = {xe C70([a,6],M") | V^x G C°([a, 6], M")}. 

In the rest of the paper, we will only consider trajectories which belong to lA. Moreover, we 
suppose that C^([a, C U, where C^([a, 6],M") is the set of C°° functions with compact 

support in [o, b]. 

First we precise the notion of dynamics in this approach. 

Definition 4. The asymmetric dynamical representation of a system S is defined by the couple 
{x+,X-) G M" X and by their respective temporal evolutions governed by the following 
differential equations 

^(x_(t),p-x^(t),...,(p-)'x_(t),t) = 0. 

The variable x+ represents the evolution in the direction past — > future, and X- in the direction 
future past. 

The time arrow is hence characterised by two objects: the global structure of the differential 
equation, via /±, and the temporal evolution operator, i.e. T>^. Whereas the direction of 
evolution clearly appears in the second one, it is not the case for the other. That is why we 
will now suppose that /+ = /_ = /. In this case, the direction of the time arrow becomes 
exclusively determined by the choice of the derivative, and the equations of the dynamics are 

f{x+{t),V+x+{t),...,{V+fx+{t),t) = 0, (2.4) 

f{x,{t),V-x.{t),...,{V-fx.{t),t) = 0. (2.5) 
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Remark 2. When = d±, the trajectory is often considered differentiahle, i.e. it verifies 
d+x{t) = d-x(t) = d/dtx{t). All information on the time arrow is therefore lost. 

In order to deal with the two directions of evolution in a unified way, we introduce the 
following functional spaces 

V = Uy.U, 
V+ = {0}, 
V" = {0} X U. 



For X = € V, we define the differential operator V by 

VX = {V^x+,V~x^). 
Consequently, for keW,V^X= {{V+fx+, {V-fx.). 

This approach will now be integrated in the framework of the embedding theories (see 

urn)- 

3. Asymmetric embedding of differential operators 

The initial motivation for the framework presented above is to conciliate fractional least 
action principle and causality. In [H HI El EO] , the Euler-Lagrange equations stemming from 
a least action principle contain both the operators and V~ . They violate the causality 
principle, as it has been noticed in [7j. For example, in [5], the equations obtained lead to 

aiL(x(t),P+x(t),t) -p-92L(x(t),P+x(t),t) = 0, 

or to 

diL{x{t),V-x{t),t) - V^d2L{x{t),V'x{t),t) = 0. 

In the framework of embedding theories, this simultaneous presence of the two operators 
is problematic for the coherence [5j: the embedding procedure and the least action principle 
are not cummutative. A solution proposed in [5] and exploited in [11] consisted in restricting 
the space of variations used in the least action principle. However, despite being very strong, 
those constraints do not however lead to a unique solution. 

A new embedding procedure, the asymmetric embedding, is presented here, and solves those 
problems. Its validity is not restricted to the fractional case. 

We begin with the differential operators introduced in [5]. 

For two vector spaces A and B, we denote T{A,B) the vector space of the functions 
f : A^ B. If f e J^(M"('=+i) X R,M'"), we define an associated operator 

F : y€.F([a,6],M"('=+i))^/(y(.),.), (g^^^ 

where /(y(»),») is defined by 

/(y(.),.) : te[a,b]^ f{y{t),t). 

If f = {fi}o<i<p and g = {gj}i<j<p are two famihes of JF(]R'^(^+1) x MjR"^), we introduce 
the operator Of defined by 



Of : xeU 



i=l 



d^ 

^(•),---,^^(-),-) , (3.2) 
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where, for two operators A = {Ai, . . . , Am) and B = {Bi, . . . , Bm), A ■ B is defined by 

(A ■ B){y) = {A^{y)Bi{y), . . . , Am{y)Bm{y)) ■ 
We extend now those operators to deal with both evolution directions. 

Definition 5. With the previous notations, the asymmetric representation of operator F, 
denoted by F , is defined by 

2 



F : {yi,y2)eJ'{[a,b] 



m{k+l) 



/(yi(*) + y2(-),-). 



(3.3) 



(3.4) 



Let Mm,2mO^) be the set of real matrices with m rows and 2m columns. We note Im the 
identity matrix of dimension m, and we introduce the operator a defined by 

X ^ {Im 0) ifXeV+\{0}, 
(0 Im) ifXGV-\{0}, 
{Im Im) otherwise. 

Now we can define the asymmetric embedding of an operator. 

Definition 6. With the previous notations, the asymmetric embedding of operator (13.21) . de- 
noted by £{0^), is defined by 

" 'F,-{V+yoG, 



£{Of) : X eV 



Fo + <^{X)Y, 



- \F, ■ {v-y o 

In particular, for (x+,0) € V+, (|3.5p becomes 



£{Of){x+,0){t) 



F^ + Y,Fi-{V^yoGi 



i=l 



and for (0, x_) G V , we have 



^(Of)(0,x_)(t) 



F^ + Y,Fi-{V~yoG, 



1=1 



X(.),...,P^X(. 



x+(t),...,(P+)VW,i 



(x^{t),...,{V-)^x.{t),t 



(3.5) 



Example 1. We set n = m = p = 1, k = 2. Let fo, fi, gi : . 
defined by 

fo{a,b,c,t) = c + e^*cos6, 
fi{a,b,c,t) = 1, 
gi{a,b,c,t) = cos a. 
The associated operator O? verifies 



be three functions 



d 



(fi ( d 

^/(^)(*) = ;772^(*) + cos -x{t) ) + ^ cos(x(t)). 



Its asymmetric embedding £{0^) is defined for trajectories (x+,x_) S V, and is given by 
£{Op{x+,X-){t) = {V+)^x+{t) + {V-)^x-{t) + e-* cos{V+x+{t) + V~x^{t)) 

'V+ cos(x+(t) + x^{t)y 



+(t(x+,x_) 



V- cos{x+{t) + X-{t)) 
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For (x-i-,0) € , the embedding becomes 

£{Of){x+,0){t) = {V+fx+{t) + e~^cos{V+x+{t))+V+cos{x+it)), 

and for (0,x_) G V~ , we have 

£{0p{0,x-){t) = {V-fx^{t)+e-^cos{V-x^{t))+V-cos{x^{t)). 

The ordinary differential equations may be written by using operators Of. Following [5], 
we consider the differential equations of the form 

Of{x) = 0, xeU. (3.6) 

Definition 7. With the previous notations, the asymmetric embedding of differential equation 
(|3.6|) is defined by 

£{Op{X) = 0, XeV. (3.7) 
Consequently, if (x+,0) S V"*", (|3.7|) becomes 



Fo + yF,-{V+yoGi 



1=1 



and plays the part of (|2.4|) . 

Similarly, for (0, x_) G V~, we obtain 



Fo + Y,F^■('^~y°Gi 

1=1 



(:E+(t),...,(P+)V(i),t) =0, 



(x_(t),...,(p-)'=x_(t),t) 



which may be related to (|2.5p . 

This method is now applied to the Lagrangian systems. 



4. Asymmetric embedding of Lagrangian systems 
We consider the same system S as above, but we suppose now that it admits a differentiable 



Lagrangian L G x M,R). In the rest of this paper, the Lagrangian L and its associated 

operator (|3.ip will be identified. For such systems, the least action principle stipulates that the 
extrema of the action provide the equation of the dynamics, called Euler-Lagrange equation. 

As we show later on, the causality problem lies in the integration by parts which appears 
in the calculus of variations. From now on, we suppose that the operators et 'D~ verify 



V+f{t)g{t)dt 



f{t)V-g{t)dt + Rab{f,g), 



(4.1) 



where f, g gU, and where Rabif,9) contains the evaluations of / and g (and possibly their 
derivatives), at points a and b. 

For example, if V+ = V- = d/dt, Rab{f,g) = f{b)g{b) - f{a)g{a). 

Concerning the reference dynamics linked to the operator d/dt, the action associated to L, 
denoted by A{L), is defined by 

A{L) : C^{[a,b],W) — > R 

X I- 



j^L(^x{t),j^xit),t^ dt, 



(4.2) 
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Its extrema provide the Euler-Lagrange equation [2] 

diL{x{t), j^x{t),t) - j^d2L{x{t), j^x{t),t) = 0. (4.3) 

Remark 3. For Lagrangian systems, the only Lagrangian determines the global structure of 
the equations. The dynamics is then completely fixed by the choice of the temporal evolution 
operator. Therefore the assumption /+ = /_ in (|2.4|) and (|2.5p seems justified for those 
systems. 

This formalism will now be fit to the approach on irreversibility by using the asymmetric 
embedding. 

4.1. Asymmetric embedding of the Euler-Lagrange equation. For X = {xi,X2),Y = 
{yi,y2) £ I^^", t € M, the asymmetric representation of the operator L, denoted by L, verifies 

L{X, Y, t) = L{xi + X2, yi + 2/2, 

A first idea consists in embedding (|4.3p directly. Given that 



dL dL N o r / \ 

- — (xi + X2,yi + y2,t) = ^ — [xi + X2,yi + y2,t) = diL{xi + X2,yi + y2,t), 
0x1 0x2 

we note diL[X, Y, t) = diL{xi + X2, yi + y2i Similarly, we note d2L{X, Y, t) = d2L{xi + 
X2,yi +y2,t). 

Theorem 1. For X £ V, the asymmetric embedding of (14. 3 p is defined by 

In particular, for (x+,0) G V^, (|4.4I) becomes 

diL{x+{t),V+x+{t),t) - V+d2L{x+{t),V+x+{t),t) = 0, (4.5) 
and for (0, x_ ) € , 

diL{x-{t),V-x^{t),t) -V-d2L{x-{t),V-X-{t),t) =0. (4.6) 

Proof. Equation (|4.3p may be written like (|3.6I) with = 1, p = 1, f = {diL, 1} and g = {d2L}. 
We conclude by using definitions [6] and [71 □ 



Those two equations may play the part of (|2.4|) and (12. 5|) for the system S. A possible 
asymmetric dynamical representation of S can hence be given by (I4.5I1 - (I4.6|) . 



4.2. Calculus of the asymmetric variations. However, it is also possible to embed the 
Lagrangian L itself instead of (|4.3p . The least action may then be applied to the embedded 
action. Indeed, the Lagrangian L may be written as (|3.2p . with p = and f = {L}. Its 
asymmetric embedding, which will be noted L, is defined by 

L : X G V I — > L{X{*),VX {•),•). (4.7) 

The associated action (|4.2p is now given by 

A{L) : V — > R 

(4.8) 



X > — > f L{X{t),VX{t),t) dt. 

J a 
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Similarly to the classical least action principle the equation of motion will be charac- 
terized through the extremal of ()4.8p . We need here to precise the notion of extremum. 

Let ^ be a vector space, B a subspace of A, and / : A ^ M a functional. Let x & A. 

Definition 8. The functional f has a B -minimum (respectively B-maximum) point at x if for 
all h € B, f{x + h) > f{x) (respectively f(x + h) < f{x)). The functional f has a B-extremum 
point at X if it has a B-minimum point or a B-maximum point at x. 

In the differentiable case, the classical necessary condition remains with this definition. 

Lemma 1. We suppose that f is differentiable. If f has a B-extremum point at x G A, then 
for all h £ B, df{x){h) = 0, where df{x) is the differential of f at x. In this case, x is called 
a B-extremal of f . 

Proof. Let h G B. The differentiable function fh ■ t i-^ f{x + th) has an extremum point in 
0. Therefore /^(O) = df{x){h) = 0. □ 

Now we introduce the space of variations 

n = G V I V/ G U, Rab{h+, f) = Rabif, /i-) = 0} . (4.9) 

Without further assumptions, we obtain the following result, derived from |51 theorem 3.11]: 

Theorem 2. Let X eV. We suppose that t ^ d-2L{X{t),VX{t),t) G U. Then we have the 
following equivalence: 

the function X gV is a TL-extremal of the action A{L) if and only if it verifies 

diL{X{t),VX{t),t) -V"d2L{X{t),VX{t),t) = 0, (4.10) 
diL{X{t),VX{t),t) -V+d2L{X{t),VX{t),t) = 0, (4.11) 

for all t G [a, b] . 

Proof. The Lagrangian L being differentiable, A{L) is also differentiable. Prom lemma [U 
X = (x+,x„) G V is a W-extremal of ^(L) if and only if for all H e H, dA{L){X){H) = 0. 
For H = (/i+, h-) G we have 

A{L){X + H) = j L{{X + H){t),V{X + H){t),t)dt, 

b 

= / L{{x+ + x- + h+ + h^){t),V+{x+ + h+){t)+V'{x- + h_){t),t)dt, 

b 

= AiL){X)+ I diLiX{t),VX{t),t){h+{t) + h^{t))dt 

b 

+ / d2L{X{t),VXit),t) {V+h+{t)+V-h^{t))dt + o{H), 

b 

= A{L){X)+ I \diL{X{t),VX{t),t)h+{t) + d2L{X{t),VX{t),t)V+h+{t)\dt 
b 

+ / \diL{X{t),VX{t),t)h-{t) + d2L{X{t),VX{t),t)V-h_{t)\ dt + o{H). 

J a 
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Consequently, 
dA{L){X){H) 



diL{X{t),VX{t),t) h+ (t) + d2L{X{t),VX{t),t) V+h+{t) 

b 



dt 



+ 



diL{X{t),VX{t),t) (t) + d2L{X{t),VX{t),t) V'h^{t) 



dt. 



Given that t ^ d2L{X{t),VX{t),t) £ U, we can use (|4.ip . No supplementary term appears 
because H £ 7i. Therefore the differential is given by 



dA{L){X){H) 



diL{X{t),VX{t),t 
b 



+ 



V'd2L{X{t),VX{t),t) h+{t)dt 
diL{X{t),VX{t),t) -V+d2L{X{t),VX{t),t)V- h^{t)dt. 



Given that X G Z^, we have t ^ diL{X{t),VX{t),t) G C°([a, 6], M*"). Consequently, 
t ^ diL{X{t),VX{t),t) -V^d2L{X{t),VX{t),t) G C°([a, 6],M'^). The space U contains 
Cg°([a, 6],M" ), so we ca n apply [ini theorem 1.2.4]: dA{L){X){H) = for all H € H if and 
only if (liTTn]l and KU\\ are satisfied. □ 

When we look at the evolution towards future, i.e. at the evolution of x+, we see that 
(x+j 0) is a 7Y-extremal of the action if and only if it verifies 

diL{x+{t),V+x+{t),t) -V-d2L{x+{t),V+x+{t),t) = 0, (4.12) 
diL{x+{t),V+x+{t),t) -V+d2L{x+{t),V+x+{t),t) = 0. 

Equation (I4.12p does not respect causality because of the operator P~. The trajectory x+ 
should hence verify two different equations, whereas only the second seems acceptable from a 
physical point of view. The same problem arises for the evolution towards past. It will now 
be shown that it is possible to overcome those difficulties by restricting the variations. 

5. Causality and coherence 

Prom a physical point of view, the elements of V are meaningless. Only their restrictions to 
V"*" and V~ are relevant. The same remark applies to Ti. Consequently, the space of variations 
has to be questioned. In (TTj, we propose to restrict the variations h by assuming V'^h = T>~h. 
Unfortunately, this hypothesis seems very strong, and may not be related to the dynamics. 
Moreover, a supplementary term appears in the Euler-Lagrange equation. In this paper we 
also restrict the variations but such problems will not arise. 

If we study the evolution towards future (we would proceed likewise for the other direction), 
it would seem natural to only consider the variations which belong also to V^. Therefore we 
introduce a new space of variations TC^ = n V^. Similarly, we set = TC Ci V~ . Erom 
theorem El we deduce the following result. 

Corollary 1. Let {x^,0) G V^. We suppose that t i-^ d2L{x^{t),'D~^x+{t),t) G U. Then we 
have the following equivalence: 

the function (x+,0) is a Ti'^ -extremal of the action A{L) if and only if x^ verifies 

diL{x+{t),V+x+{t),t) - V-d2L{x+{t),V+x+{t),t) =0, (5.1) 

for all t G [a, b] . 
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We obtain a single equation, which is precisely (|4.12|) . the problematic one. Consequently, 
the variations (in Ti) and the trajectories (in V) cannot have the same status. 
On the contrary, if the trajectories are chosen in H~ , the problem is solved. 

Theorem 3. Let G V+. We suppose that t ^ d2L{x+{t),'D+x+{t),t) G U. Then we 

have the following equivalence: 

the function (x+,0) is a 7i~ -extremal of the action A{L) if and only if x^ verifies 

diL{x+{t),V+x+{t),t) - P+92L(x+(t), P+x+(t), t) = 0, (5.2) 

for all t £ [a, b] . 

We obtain once again a unique equation, but which is now causal in the direction past — > 
future. We have a similar result for the other direction. 

Corollary 2. Let (0, x_) G V~. We suppose that t i-^ d2L{x^{t),'D^X-{t),t) G U. Then we 
have the following equivalence: 

the function (0, x_) is a Ti.'^ -extremal of the action A{L) if and only if x^ verifies 

diL{x-{t),V- x-{t),t) - V- d2L{x-{t),V- x-{t),t) =0, (5.3) 

for all t G [a, b] . 

This equation is causal in the direction future past. 

Both of the equations (|5.2I) and (|5.3p stem from a least action principle, and moreover 
respect the causality principle. In addition, no significant restriction is done on the non-zero 
components of the variations. Therefore we obtain a least action principle similar to the 
classical one (with d/dt), except that trajectories and variations are not ruled by the same 
dynamics. 

This result may seem surprising: it shows that the equations of the dynamics in a given 
temporal direction are obtained through variations evolving in the opposite way. Let us discuss 
this paradox. 

The least action principle is a global vision of the dynamics: the trajectory is directly 
determined on its whole temporal interval [a, 6]. According to Poincare, the system "seems 
to know the point to which we want to take it". In [13], an history of this principle is 
presented, and it is shown that this finalist aspect has been the subject of controversies since 
its formulation by Maupertuis in 1746 [H]. By making the past and the future of the trajectory 
appear explicitly, equation (|5.ip is in agreement with this global approach. For example, in the 
fractional framwork, V'^ = aX'f and takes into account the whole past ([a,t]) of the trajectory 
and P~ = — t'L'^ , all the future ([t, ?)]). So it may seem that variational formulation (global) 
and causal equation could be incompatible. But as it has just been shown, the variations can 
lift this difficulty. Because they obey to the reverse dynamics, they "catch" the anti-causal 
part of the least action principle. Then the equations of the dynamics can respect causality. 
In the sum x+ + we add two comparable functions, but their underlying dynamical and 
physical natures differ. The trajectory may be called "real", actual, while the variation h- 
may be seen as "virtual", potential. To sum up, the finalist aspect of the least action principle 
may lie in the nature of the variations. Given that these do not possess a concrete realisation, 
this problematic characteristic of the principle may seem less disconcerting. 

Remark 4. In the case of a dynamics governed by d/dt, this discussion is obscured because 
of the local aspect of the derivative. 
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Furthermore, we note that (|5.2p and (|5.3|) are identical to (|4.5|1 and (|4.6I) . So we obtain the 
following commutative diagram 

L{x{t), ix{t),t) — L{X{t),VXit),t) 



LAP 



LAP 



[diL - id2L] ix{t), ix{t),t) = — [diL - V^d2L] {x±{t),V^x±{t),t) = 0, 

where LAP means "least action principle" and AE "asymmetric embedding". Following [5], 
we say that the asymmetric embedding is coherent. 

This "robustness" motivates the adaptation of definition H] for Lagrangian systems. 

Definition 9. The asymmetric dynamical representation of a Lagrangian system with La- 
grangian L is defined by the couple G M"' x and by their respective temporal 
evolutions governed by the following differential equations 

diL{x+{t),V+x+{t),t) -V^d2L{x+{t),V^x+{t),t) = 0, 
diL{x-{t),V'x-{t),t) -V~d2L{x^{t),V'x^{t),t) = 0. 

The asymmetric embedding of those systems is now illustrated with few operators V'^ and 

V-. 



6. Particular cases 

We consider the same Lagrangian system <S as above, with Lagrangian L. First the case 
of the finite differences is presented. These provide an clear illustration of the topic. Then 
the degenerate case of the classical derivative is adressed, before finishing with the fractional 
operators. 

6.1. Finite diff'erences. For e > fixed, we choose I?^ = df, with 

dtf(t) = (6.1) 
djfit) = -^{f{t + e)-f{t)). 

We verify that the operator df takes into account the past of / and d^ the future. 
For a = -oo and 6 = +oo, we have U = CO(M,M"). We verify that C;?°(M,M") C U. 
We also have a relation similar to (|4.ip : 

/ dtf{t)g{t)dt = - [ f{t)d-g{t)dt. 
In this case, Rabif, ff) = and H = V. 

For those operators, the asymmetric dynamical representation is 

diL{x+ (t) , dfx+ (t) , t) - dfd2L{x+ (t) , dtx+ (t) , t) 
diL{x^{t), d^x^{t),t) — d~d2L{x^{t),d~x^{t),t) 



0, 
= 0. 



(6.2) 
(6.3) 



Remark 5. Concerning the least action principle, we may note that the condition t i— > 
d2L{X (t) ,VX (t) , t) is always verified in this case. 
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This example is a good illustration of the embedding notion. When we observe (or simulate) 
the evolution of a classical Lagrangian system (ruled by the operator d/dt), we only measure its 
state at particular instants. The observations remain ponctual, never continuous. The velocity 
is calculated from those ponctual measures, with a formula similar to (|6.1|) . Experimentally, 
we cannot choose e as small as we want in order to recover d/dt. Besides, even if we suppose 
the time step very small, we would approximate d+ defined by (|2.2p and not d/dt. The 
dynamics that we observe is therefore not ruled by d/dt, but by d^. 

6.2. Classical derivative. Nevertheless, let us look at the case where the embedding does 
not modify the temporal evolution operator. We have = D" = d/dt. We consider an 
interval [a, with -co < a < b < +00. In this case, U = C^{[a,b],W^). Once again, 
C^(M,M") C U. The relation (|4.ip is the classical integration by parts 



Here we have Rab{f^g) = fib)g{b) — f{a)g{a) and Ti = {(/i+,/i_) € V|/i+(a) = h^{a) = 



The equations are identical in both directions, are moreover similar to the reference equation 
(14. 3p . According to our approach on irreversibility, the system is in this case reversible. This 
point could be discussed. 

6.3. Fractional derivatives. For a detailed presentation of the fractional calculus, we refer 
to [TSl [m El] . For fractional derivatives, the role of the past and the future clearly appear. 
Regarding the different non causal Euler-Lagrange equations derived in [H HI O [20] , the main 
contribution of this paper is the obtention of a causal equation. 

We consider an interval [a, 6], with —00 < a < b < +00, and we consider the Riemann- 
Liouville fractional derivative of order a, with < a < 1. For the fractional case, precising U 
is not easy. We refer to [21] for precisions. From [2TJ p. 159], we still have C^(M,M") C U. 

Moreover, as indicated in [12], we introduce a extrinsic time constant r in order to retain 
the dimensional homogeneity of the equations. By the way, let us mention that the methods 
exposed in [12] naturally apply to the asymmetric embedding. Therefore it is possible to 
obtain fractional equations which stem from a least action principle, and which additionnally 
preserve causality and homogeneity. So we set = r"~^ and T>~ = — r"~^ , with 



Those operators deal with the integrality of the past and the future of the function. The 
question of causality is hence momentous here. Relation (14. ip is now 




h+{b) = h-{b) = 0}. 
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r / \ '2 2 2 
L{X, V) = - V — CO X . 



which imphes Rabif,g) = and Ti = V. If we had chosen the Caputo fractional derivative, 
the only difference would have been Rab{f,g) / 0. 

The asymmetric dynamical representation is given here by 

5iL(x+(t),r"-\Pf -r°-\Pf a2L(x+(t),r°-\Pf = 0, (6.4) 

aiL(x_(t), -r"-i x.{t),t) + r"-i tV^ d2L{x^{t), -r^-^ x.{t),t) = 0. (6.5) 

Example 2. We finish with the case of the harmonic oscillator, with a = — oo and b = +oo. 
Its Lagrangian may be written as 

1 
2 

In this case, (16. 4p and (16. 5|) are 

T^^-^{^^Vf o ^^Vf)x+{t)+uj^x+{t) = 0, 
T^^-\tV%^o,Vl^)x^{t)+u^x.{t) = 0. 
Given that -oo^^f o _^Vf = _oo^??" and tT^X^ o tV^^ = t'D'^oo , we obtain 

T^''-^-oo'D^'^X+{t)+LO^X+it) = 0, 

r2"-^p2a^x_(t)+c^2^_(t) = 0. 
For a = I, _ooI^t" = tT^+^oo — /dt^ . We recover the classical equations 

-^x+{t)+uj^x+{t) = 0, 

■^x-{t) +uj^x^{t) = 0. 

The dynamics is reversible, which is in agreement with the usual notion of reversibility 
(invariance under the transformation t^ —t). 

On the other hand, for a = 1/2, we have ^^cV't" = d/dt and t^'^ao — —d/dt, which leads 

to 

j^X+{t)+TUj'^X+{t) = 0, 
-^X^{t) +TU?X-{t) = 0. 

Those two equations differ: the dynamics is therefore irreversible. Once again, this result 
agrees with the usual definition of irreversibility. 

7. Conclusion 

By distinguishing left and right derivatives, a formal approach on irreversibility has been 
presented. It may help clarify the role of past and future in the least action principle. In 
particular, trajectories and variations have been uncoupled at the dynamical level. This 
differentiation has led to the obtention of Euler-Lagrange equations which respect causality, 
even in the fractional case. 

We underline once again that this point of view on irreversibility is formal and does not 
explain the physical origin of this phenomenon. However, it may constitute a trail for a better 
understanding, through the role of the temporal evolution operators. Without closing the 
debate on the least action pinciple, this paper precises the relation between this principle and 
causality. Moreover, the formalism presented here associated to that in [T2] finally conciliates 
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fractional Euler-Lagrange equations with fundamental physical principles. We hope that this 
mathematical framework, now physically satisfying, will facilitate the description of phenom- 
ena still not properly understood, particularly in the fractional case. We think for example 
of some Hamiltonian chaotic systems, for which fractional diffusion models have been devel- 
opped in [25 . In those systems, fractional dynamics might arise with the construction of a 
new "macroscopic" time |S], based on recurrence times and describing long-term behaviours. 
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